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Thin discs, thick dwarfs and the effects of stellar feedback 
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ABSTRACT 

We investigate the role of stellar mass in shaping the intrinsic thickness of galaxy 
discs by determining the probability distribution of apparent axis ratios (b/a) for two 
different samples that probe the faint end of the galaxy luminosity function. We find 
that the b/a distribution has a characteristic 'U-shape' and identify a limiting mass 
M* w2x 10 9 M below which low-mass galaxies start to be systematically thicker. 
This tendency holds for very faint (Mb ~ —8) dwarfs in the Local Volume, which are 
essentially spheroidal systems. We argue that galaxy shape is the result of the com- 
plex interplay between mass, specific angular momentum and stellar feedback effects. 
Thus, the increasing importance of turbulent motions in lower mass galaxies leads 
to the formation of thicker systems, a result supported by the latest hydrodynamical 
simulations of dwarf galaxy formation and other theoretical expectations. We discuss 
several implications of this finding, including the formation of bars in faint galaxies, 
the deprojection of Hi line profiles and simulations of environmental effects on the 
dwarf galaxy population. 

Key words: galaxies: formation - galaxies: structure - galaxies: spiral - galaxies: 
dwarf - galaxies: statistics 



1 INTRODUCTION 

Galaxy discs are key structural components in the under- 
standing of galaxy formation. They con tain approximately 
60% of the stellar mass in the Universe IjDriver et alj|2007h 
and ar e one of the ma in sites of current star formation ac- 
tivity (|Kennicuttl 199a ) . Moreover, their prominence and ap- 
pearance form the basis of the Hubble ( f 926) morphological 
sequence. 

Despite this enormous importance, understanding the 
details of disc formation remains extremely challenging. 
The basic scenario considers that baryons cool and collapse 
within hierarchically assembled dark matter haloes. A cer- 
tain amount of their torque-acquired angular momentum is 
transferred to the baryonic component lea ding to the forma- 
tion of a rotationally supported thin disc (jFa.ll fc Efstathioul 
If980l ). This relatively simple picture has been able to re- 
produce many observational properties of disc galaxies - 
flat rotation curves, Tully-Fisher relat i on, gas content (e.g. , 
Dalcanton et al.llf997l : iMo et al.lli~99Sl : Ivan den Bosch] If 9981 . 



200(3 )- but suffers from several shortcomings. On the one 



hand, detailed N-body simulations showed that when the 
dissipative effects of gas are not considered, the hierar- 
chical nature of structure formation unavoidably results 
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in disc destruction l|Toth fc Ostrikerllf992l ). On the other 
hand, hydrodynamical simulations of this process system- 
atically produced discs that were too small and too cen- 
trally concentrated due to excessive angular momentum 
exchange between the gas and the dark matter haloes 
(JNavarro fc SteinmetJI 19971 ) . 

A physically-motivated solution for these problems 
is the inclusion of strong feedback effects from differ- 
ent sources -but specially from star formation and super- 
novae explosions- which, when coupled with a cosmic UV 
field, are able not o nly to produce realistic discs (e.g., 
iGovernato et al.ll20fd ). but also to provide on e possible ex- 
planation for the missing satellite problem (jKlypin et al.l 
1 19991 ) . 

This approach, in turn, raises the question of which 
role does mass play in shaping the properties of discs, as 
these strong heating mechanisms are exp ected to produce 
a gre ater influence in lower-mass galaxies (jKaufmann et al.l 
120071 . KWB07 hereafter). Indeed, it is well known that galax- 
ies h ave higher gas mass fractions (e.g., ISchombert et al.l 
2001) and more extende d star formation time scales 
( Hunter fc Gallagher! If 98a ) as they are less massive, and it 
has been proposed that dwarf an d disc galaxies ar e probably 
two different structural entities (|Schombertll2006l) . 

From an observational point of view, galaxy discs are 
best described as flattened triaxial ellipsoids with exponen- 
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tial surface brightness profiles (|FreemarJI 19701 ) . Though they 
are traditionally considered to be perfectly circular, it is well 
known that dis cs are indeed slightly elliptical, with b/a > 0.9 
(jLambas et al.lll992l ; lRvdenll2004r ). They are said to be thin 
because their vertical to radial axis ratios have long been 
known to lie in a narrow r ange 0.15 < go 5, 0-25 IjHolmberg 
Il95d ; ISandaae et aflll970l : iHeidmann et alJll972h . However, 
this thickening is not unifor m, as late-type spirals have thin- 
ner discs than early-types (jBottinelli et al.l 1 19831 ; I Guthrie! 
1 19921 ). If we neglect these small ellipticity deviations, discs 
can then be thought of as axisymmetric oblate spheroids 
with intrinsic thickness go. In that case, the apparent axis 
ratio of a disc is determined by b/a = i/g 2 , + (1 — go) 2 cos 2 i, 
where i is the corresponding inclination angle. Interestingly, 
the apparent axis ratio distribution of a randomly oriented 
population of such discs peaks at a minimum value of go _ 
i.e., edge-on galaxies are much more common than any other 
projection. This property implies that the intrinsic thickness 
of a disc population can be identified as the minimum value 
of the distribution of apparent axis ratios. 

In this paper we take advantage of this prediction to 
investigate the range of masses where thin discs exist. For 
this purpose we study the probability distribution of b/a 
as a function of galaxy mass (luminosity) for two different 
samples that probe the faint end of the galaxy luminos- 
ity function. In Section 2 we develop further on the sam- 
ples characteristics and the methodology. Section 3 presents 
the main results, which are discussed in Section 4 together 
with additional implications of this investigation. Through- 
out this Letter we adopt a ACDM cosmology with flo = 0.3, 
fl A = 0.7 and h n = H /(75 kms -1 Mpc -1 ). 



2 SAMPLE AND METHODOLOGY 

Two different galaxy samples have been used throughout 
this work. The first sample consists o f all 9245 galaxies in 
the SDSS-DR7 (jAbazaiian et alj|2009j ) with recession veloc- 
ities in the 2000 < cz < 6000 kms -1 range. The lower limit 
was set to avoid strong correcti ons in galaxy distances due t o 
the Virgo infall velocity field (jTammann fc SandageHl985T ). 
The upper one ensures that the sample is volume-limited 
for galaxies brighter than Mi ~ —16.5, additionally includ- 
ing 2063 fainter objects down to a limiting magnitude of 
Mi w —14.5. We however note that this incompleteness only 
affects the statistical significance, but not the shape, of the 
axis ratio probability distribution at the lowest luminosities, 
as the volume we probe is too narrow (25 < D < 80 Mpc) for 
any shape evolution to take place. A further concern might 
arise given the well-known incompleteness of t he SPSS imag- 
ing fo r low surface brightness galaxies (LSB; iBlanton et al.l 
2005). However, if a population of LSB discs exists, we will 
preferentially miss the face-on (high b/a) systems instead of 
the edge-on, higher surface brightness discs. Therefore, the 
fraction of highly flattened galaxies is in any case a lower 
limit. 

The SDSS imaging provides ugriz magnitudes and sev- 
eral estimates of axis ratios in each band. For this study 
we have used the 25 magarcsec -2 isophotal axis ratios, 
which provide a robust measurement of galaxy shapes at 
their outer regions -a few times the galaxy's effective radius 
|Vincent fc Rvdenll2005l )-. where the effects of dust and the 



presence of bulges are less important. Seeing effects do not 
play a role either, as the proximity of the sample ensures that 
essentially all galaxies have minor axis sizes larger than the 
typical seeing FWHM, with a median of 20 arcsec in the 
r-band. SDSS isophotal b/a can sometimes take unrealistic 
values in a given photometric band. To avoid these, we have 
used the median of the three high signal-to-noise gri axis ra- 
tio estimates as our reference b/a value, therefore producing 
a less noisy (o b / a « 0.05) shape estimate for each galaxy. 
We however verified that our results are independent of the 
selected filter. 

Given that the main purpose of this work is to study the 
influenc e of galaxy mass in the formation of discs, we have 
used the lBell et all IJ2003h M/L ratios with a lKroupal (|200lh 
initial mass function in the transformation from luminosities 
to stellar masses. 

In order to extend our analysis to even fainter luminosi- 
ti es, we also analyse the cat alogue of neighbouring galaxies 
of lKarachentsev et al.l 1)20041 ). This is an all-sky catalogue in- 
cluding 445 galaxies brighter than Mb < —8 with individual 
distance estimates D ^ 10 Mpc and is expected to be 70%- 
80% complete within 8 Mpc. From this catalogue we used 
the B-band absolute magnitudes and the 25 magarcsec -2 
isophotal apparent axis ratios b/a of each individual galaxy. 



3 RESULTS 

In the left panel of Fig.[T]we present the probability distribu- 
tion (grey scale) of apparent axis ratios in intervals of stellar 
mass for the SDSS sample. This figure contains several in- 
teresting features. 

First, the fact that, aside from the most massive bins, 
the 99% percentile (upper dot-dashed line) slightly deviates 
from unit y, consistent with the small ellipticity typical of 
discs f e.g. . lRvdenl2004r ) . Second, and most important for our 
purposes, both the median and the 1% percentile (solid and 
lower dot-dashed line, respectively) of the b/a distribution 
present a characteristic 'U-shape\ with a minimum b/a w 
0.15 located at M* fa 2 x 10 9 M Q B We note that this value 
corresponds to Mi ~ — 18, a much brighter luminosity than 
the limiting magnitude of our volume-limited sample. The 
minimum in the b/a distribution is better appreciated in 
Fig. [T] (middle), where we plot the fraction of galaxies having 
an axis ratio smaller than a given value as a function of 
stellar mass. The distribution is indeed a strong function of 
M* , and the fraction of low- mass thin galaxies never exceeds 
a few per cent. 

The increase of apparent axis ratio in the high mass 
range can be ascribed to two effects. The main reason is that 
the contribution of the spheroidal component is more impor- 
tant at higher M« , being totally domi nant for the brightest 
ellipticals fcf. Ivan der Wei et al.l 120091 ). Additionally, early- 
type (more massive) spirals tend to have thicker discs than 
their late-type counterparts (Bottinclli ct al. 1983). 

The increasing thickness towards lower masses is per- 
haps more puzzling given the previous fact and the decreas- 
ing bulge mass fraction at faint luminosities. Galaxies with 



1 Recall that, as previously pointed out, the lower percentile can 
be interpreted as the intrinsic thickness of the flattest galaxy pop- 
ulation in a given stellar mass interval. 
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Figure 1. Left: the grey scale represents the probability distribution of apparent axis ratios in intervals of stellar mass for the SDSS 
sample -where dark (light) indicates a high (low) probability within each M* bin. The median (solid line) and the 1% and 99% percentiles 
(dot-dashed lines) are also indicated. The distribution has a characteristic 'U-shapc', showing a trend of increasing axis ratios for lower 
(and higher) mass galaxies. Recall that the typical b/a uncer t ainty is RiO.05, i.e., the binsize. The starred symbols indicate the edge-on 
thickness of the dwarf galaxies simulated bv lGovernato et al.l 42010C) —and are therefore lower limits for their observed b/a. Middle: each 
curve shows the fraction of galaxies in each mass interval having a lower b/a than the specified value. The peak of the distribution occurs 
at a characteristic mass M* ss 2 X 10 9 Mq, which we identify as the minimum stellar mass of thin disc galaxies. Right: same as in the 
left panel but for the Local Volume sample, that includes galaxies down to Mg = —8. Note the lack of faint flattened systems. 



stellar masses M* ~ 10 Mq are most probably spheroidal 
systems, with apparent axis ratios in the 0.3 < b/a < 0.95 
range and a median b/a ~ 0.8. This paucity of thin discs 
in the dwarf regime is enhanced at the faintest magnitudes 
probed by the Local Volume catalogue (Fig.Q] right), with 
less than 1% of the Mb > — 14 galaxies having b/a < 0.3. 



4 DISCUSSION AND CONCLUSIONS 

Hints of incre asing thickness for less massive galaxies are 
not n ew (e.g.. lrleidmann et alj|l972t lYoachim fc Dalcantonl 
120061 ) , but to our best knowledge we are the first to identify 
a characteristic mass below which low-mass galaxies start to 
be systematically thicker. The question then arises naturally. 
What causes this effect? 

One possibility is that environmental effects could in- 
fluence galaxy shape. While the narrow volume we probe 
does not include any significant massive cluster -it is just 
further than Virgo and closer than Coma-, there remains 
the possibility that low-mass galaxies are tidally affected 
by more massive companions in group environm ents, as has 
been previously observ ed in the Local Group |Choi et al.l 
120021 ; iGeha et all I2OO6I ) . To test for this scenario, for each 
SDSS galaxy we h ave computed a tidal parameter (e.g.. 



xy 

a [2 



IVarela et all I2004T I T P = max[(M», P /M»,0 (A 25 /D P ) 3 ], 
where M»,j and A25 are its stellar mass and isophotal ma- 
jor axis, and M*,p and Dp are the corresponding stellar 
mass and projected distance to any potential perturber hav- 
ing a relative velocity |Av| < 1000 kms -1 and a lumi- 
nosity Mi < —16.5. This parameter is proportional to the 
maximum ratio of external to internal forces that act on a 
galaxy, and takes into account that the effects of a closer, 
less massive companion can be even more important than 
those originated by a massive, more distant galaxy. Fig. [2] 
shows in grey scale the probability distribution of Tp in 
each mass interval for the SDSS sample, while the different 



lines indicate the corresponding 1%, 10%, 50%, 90% and 
99% percentiles. There is a clear trend of decreasing Tp 
with decreasing mass, i.e., low-mass galaxies are in general 
less tidally affected than their more massive counterparts. 
This result is consistent with the fact that the galaxy lumi- 
nosity fu nction is dominated at all luminosities by central 
galaxies (|Coorav fc Milosavlievicll2005l ). Given that the nor- 
malisation of the tidal parameter is somewhat arbitrary , we 
have used the Hyperleda database l|Paturel et al.ll2003l ) to 
compute the Tp of three dwarf galaxies experiencing differ- 
ent l evels of interaction with M31 (see iKarachentsev et al.l 
[200J). We find that the fraction of strongly tidally perturbed 
galaxies like M32 or N205 is only ~10% , while the vast ma- 
jority of low-mass galaxies are rather isolated systems as 
WLM. 

We therefore conclude that tidal effects are not respon- 
sible of the increasing thickness found in fainter galaxies, and 
suggest that the effect is probably related to the increasing 
importance of feedback mechanisms in low-mass haloes. 

KWB07 investigated the effect that a gas temperature 
floor (Tp) -as might arise in the presence of a cosmic UV 
field and/or due to stellar feedback- has on the properties 
of low-mass galaxies. Their instructive fig. 1 shows that the 
final morphology of a galaxy is the result of a complex in- 
terplay between the turbulent support provided by this Tp 
and the angular momentum support, characterised by the 
dimensionless spin parameter A. For low-mass haloes, the 
pressure support radius becomes comparable to or larger 
than the rotational support radius, and thus galaxies are 
naturally formed thicker. 

ISchomberj (|2006l ) was the first to note a separation 
between dwarf and disc galaxies in the scale-length versus 
mass plane, in the sense that the former are more diffuse 
(extended) than the latter. He suggests that the increasing 
importance of turbulent motions in dwarf galaxies results in 
thicker systems, potentially giving rise to the two observed 
sequences. However, he also points out that 'this scenario 
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Figure 2. Probability distribution of the tidal parameter in bins 
of stellar mass (grey scale) for the SDSS sample. Note that low- 
mass galaxies are in general less tidally affected than their more 
massive counterparts. For reference we show the corresponding 
values for three local- volume dwarf galaxies experiencing different 
levels of interaction with M31. 



provides no explanation for why there are two distinct se- 
quences rather than a smooth transition in structure from 
an oblate shape to a triaxial one as one progresses to lower 
galaxy mass'. We believe that the key to this separation re- 
sides again in the complex interplay between galaxy mass, 
heating mechanisms and specific angular momentum. 

Going back to KWB07's fig. 1, we can identify three dif- 
ferent regions. For the most massive haloes, angular momen- 
tum support prevails regardless of spin parameter or temper- 
ature floor value. In this region, thin discs will form, lying 
in Schombert's sequence defined by Sc-Sd galaxies. In the 
case of very low-mass haloes, the opposite is true: pressure 
support is always dominant even for haloes with the highest 
spins A ~ 0.1. Galaxies in this region will be puffier, creat- 
ing the Im and dlrr sequence in the scale-length vs. mass 
relation. Finally, for the intermediate mass region KWB07 
point out that a range of morphologies is predicted if galax- 
ies have spins with in the expected range (0.01 < A < 0.1; 
iBullock et al.ll2001D and the temperature floor does not vary 
too much among galaxies. This would naturally explain the 
co-existence of the two sequences in the narrow 10 9 -10 10 M© 
mass range: galaxies in the high-spin tail of the distribution 
would end up in the disc sequence, while the dwarf sequence 
would be populated by galaxies with low spin parameter. 

Evidences for the decreasing importance of ordered mo- 
tions in favour of turbulent ones come from the high thick- 
ness of the neutral gas co mponent in very faint gala xies 
from the FIGGS sample (|Rovchowdhurv et al.1 12010| '). or 
fr om the shape of H i line profiles from the dwarf sample 
of iGeha et al.l (|2006l ) . where onl y 18% of the obje cts exhibit 
double-horned profiles (but see ISchombertJ I2006T ) . Further- 
more, this fraction increases to 30% if only edge-on galax- 
ies are considered, suggesting that the shape of these flat 
objects can be driven by rotation. Indeed, the projected el- 
lipticity of a rotationally supported oblate spheroid with a 
constant anisotropy parameter is an increasing function of 



v/a i|Binnevlll978r ). so that objects with prevailing ordered 
motions are expected to be more flatte ned. This picture is 
supported by fig. 8 of IGeha et al.l (2006), where it is shown 
that the small number of dwarf galaxies with low apparent 
axis rati os always lie above the best-fit baryonic Tully-Fisher 
relation IjMcGaugh et al.ll2000l ). i.e., their rotation is higher 
than the mean for their baryonic masses. 

How do all these results compare with ours? The mini- 
mum of our b/a distribution occurs at a characteristic mass 
(luminosity) M* « 2 x 10 9 M (M % ~ -18). This is roughly 
the same value as the lower boundary of Schombert's disc 
sequence, and would thus indicate a limiting mass for thin 
discs. The baryonic Tully-Fisher relation derived by Geha 
et al. indicates that galaxies below these luminosities reside 
in haloes with V c < 80 km s _1 and this suggests that if their 
thickness is set by turbulent motions, the associated tem- 
perature floor must be high (see fig. 1 of KWB07) . 

In order to check if stellar feedback effects can be ef- 
fectively responsible of the previous trend, we have derived 
the intr insic thickness o f the d warf galaxies recently mod- 
elled bv lGovernato et al.l (120100 . These hydrodynamical sim- 
ulations have first produced realistic dwarf galaxies thanks 
to their high resolution and a detailed treatment of bary- 
onic processes. Their simulated DG1 has a stellar mass 
M* = 4.8 x 10 8 Mq and a dimensionless spin parameter 
A = 0.05. We have measured a 25 magarcsec -2 isophotal 
axis ratio of qo — 0.3 using an r-band edge-on image of 
the dwarf (five-pointed star in Fig. 1, left). This value is 
higher than the minimum thickness we obtain for similar 
mass galaxies, and we speculate as to whether these flatter 
objects could reside in haloes with higher spin parameters. 
Interestingly, the second galaxy simulated by Governato et 
al. (DG2) -which is slightly less massive, has the same A and 
had a quieter merging history than DG1-, is even thicker 
(qo — 0.35, see Fig. 1), nicely following the trend with mass 
we find. We therefore suggest that star formation feedback 
effects -which can remove low angular momentum material 
and produce bulgeless dwarfs with shallow central dark mat- 
ter profiles- are also responsible of the increasing thickness 
of low mass galaxies. The shallower potential well certainly 
allows more turbulent (and therefore also vertical) motions, 
so stars naturally settle in a thicker disc. This scenario is 
additionally appealing as it can provide an explanation (cf. 
KWB07) for the observed high gas mass frac tions and low 
star formation efficiencie s of faint systems (jWarren et al.l 
I2007J : lAmon'n et alj|2009h . 

The results presented in this study have several further 
implications: 

1) The existence of a minimum stellar mass for thin discs 
has a pronounced impact on the formation of bars, as 
these structures are thought to be formed through instabil- 
ities in cold (thin), rotationally supported discs. Therefore, 
bar formation could be inhibited in galaxies less massive 
than M* « 2 x 10 9 Mq simply because they are not cold 
en ough. This is consistent with the result recently obtained 
by iMendez-Abreu et al.l (J2010j), which show that bars are 
only found in the narrow 10 9 -10 11 M q mass range. 

2) It has long been noticed (e.g., iBottinelli et al.l [1983) 
that a varying intrinsic thickness can introduce serious sys- 
tematic errors in calculated inclinations if they are com- 
puted assuming discs have a universal, fixed go -as it is 
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actually usually done when deprojecting Hi line profiles. 
While errors are only of a few per cent if the galaxy is 
still an oblate spheroid, the situation becomes more criti- 
cal if dwarf s are actually bet ter represented as triaxial ellip- 
so ids (e.g., Sung et al.lll99a). Moreover, the recent results 
of iRovchowdhurv et al.l (j2010f l indicate that gas discs are 
also thicker in faint galaxies, further increasing the previous 
concerns. 

3) As already pointed out by KWB07, simulations that 
investigate the role of environmental effects on the popula- 
tion of dwarf galaxies usually consider a low-mass thin disc 
as their fiducial model. Given that real dwarf galaxies are 
intrinsically thicker -and have shallow central dark matter 
profiles-, they are likely more susceptible to dynamical in- 
teractions than most simulations predict. 

Interestingly, the evidence of increasing thickness with 
decreasing mass can be found in several other sam- 
ples of field late-type dwarfs in the literature (e.g., 
lEder fe Schombertl l200d : iGeha et all l2006h as well as fo r 
non-nucleated dEs in the Virgo cluster IjLisker et al.l l2007). 
further suggesting that it might be related to a formation 
process independent of environment. We believe that future 
and detailed hydrodynamical simulations of dwarf galaxy 
formation -probing a variety of masses, spin parameters and 
merging histories- should be able to reproduce the proba- 
bility distribution of apparent axis ratios and provide more 
quantitative explanations of the observed trends. 
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